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Abstract. The dynamics of metric perturbations is explored in the grav- 
ity theory with anomaly-induced quantum corrections. Our first purpose is 
CJ . to derive the equation for gravitational waves in this theory on the general 

^ ■ homogeneous and isotropic background, and then verify the stability of such 

^ background with respect to metric perturbations. The problem under consid- 

eration has several interesting applications. Our first purpose is to explore the 
stability of the classical cosmological solutions in the theory with quantum 
effects taken into account. There is an interesting literature about stability 
of Minkowski and de Sitter spaces and here we extend the consideration also 
to the radiation and matter dominated cosmologies. Furthermore, we analyze 
the behavior of metric perturbations during inflationary period, in the stable 
phase of the Modified Starobinsky inflation. 
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1 Introduction 



The semiclassical approach to gravity is usually associated with the equation 

Gfj,u = Rfiu — 2 R9^lv = < T^i, > (1) 

and implies that the gravity itself is not quantized. The averaging in the r.h.s. of the 
last equation comes from the quantum matter fields. Different from what one can think, 
the r.h.s. may be very nontrivial even if no matter sources are present in the given point 
of space-time. The reason is that the average of the Energy-Momentum tensor of the 
vacuum can be nontrivial, for it may depend on the curvature tensor components and 
its derivatives, with possible non-local structures. There are two subtle points in the 
equation ([1]). Let us start from the terminology. The renormalizable theory of matter 
fields on curved space-time background requires that the action of gravity should be 
extended compared to the one of General Relativity (GR) [H E], [3] (see also |1] for the 
recent review). The full action includes Einstein-Hilbert term, which is the origin of the 
r.h.s. of ([1]) with the cosmological constant term 

Seh = - J d'x^ {R + 2K) (2) 

and also the higher derivative terms 

Shd = j d^xy/^ {aiC^ + a2E + a3DR + aiR'^} , (3) 

where = Rl^^i^ ~ '^Ra/s + (1/3) R^ is the square of the Weyl tensor and E = -R^^„^ - 
4i?^^ -|- R^ is the integrand of the Gauss-Bonnet topological term. All terms of the action 
of vacuum 

Svac = Seh + Shd , (4) 

belong to the gravitational action. If we consider the Einstein equations as something 
intended to define the relation between geometry and distribution of matter, it is clear 
that the whole Syac should contribute to the l.h.s. of (JT]). However, by traditional virtue 
we use to put all the contributions beyond the Einstein tensor to the r.h.s.. To a great 
extent this way of settling the terms is explained by the fact that Eq. ([1]) with the 
classical Energy- Momentum tensor at the r.h.s. work pretty well and provide very 
good fit for many observational tests of GR. So, it looks like all extra terms, except the 
cosmological constant, are in fact unnecessary at the classical level and therefore their 
introduction is no justified in a purely gravitational framework. In other words, why 
should the gravitational physicist who works with the very large scale phenomena, worry 
about the quantum notions, such as renormalizability? In fact, the use of the terms 
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([3]) may lead to serious problems, because these terms are known to produce unphysical 
ghost terms for the linearized gravitational field on the fiat background [5]. So, it is 
somehow unclear how to apply the consistent at quantum level theory (j4j) for the classical 
gravitational purposes. 

Another aspect of the same problem is related to quantum corrections to (j4]). In gen- 
eral, the problem of deriving such corrections is unsolved, but there is one important case 
where the situation is quite clear. The one-loop effective action of massless conformal 
fields is essentially controlled by conformal anomaly [6]. Indeed, the anomaly- induced 
effective action of vacuum includes an arbitrary conformal functional, but in many par- 
ticular cases its role is known to be very restricted. For instance, this functional is a 
trivial constant for the homogeneous and isotropic metric, where the anomaly-induced 
quantum corrections produce Starobinsky infiation [7]. In the case of black holes, taking 
into account the conformal anomaly enables one to calculate Hawking radiation |8] and 
using the anomaly-induced action one can even classify the vacuum states in the vicinity 
of the black hole P, [TU] . The equations for gravitational waves calculated by using direct 
methods [HI |12] and anomaly-induced action [13\ produce equivalent results at least on 
the de Sitter background. In both last cases the mentioned conformal functional appar- 
ently plays no role, which can be explained by the fact that the conformal anomaly picks 
up all quantum effects which correspond to the UV limit and, therefore, the remaining 
terms can be relevant only for the sub-leading effects. So, from the Quantum Field The- 
ory viewpoint the anomaly-induced effective action of vacuum represents a well-defined 
quantum contribution which can be used to verify the compatibility of gravity and quan- 
tum effects in the sense we have discussed above. The late Universe represents, in fact, 
a very good opportunity for such a verification. The typical energy scale of gravity is 
given by the Hubble parameter or by the scale of the gravitational waves, both of which 
are much smaller than the masses of all quantum massive fields. Certainly, all such fields 
strongly decouple in the higher derivative sector of the theory [H] and therefore the only 
active quantum field is photon. So, we arrive at the conclusion that the anomaly- induced 
action of vacuum coming from the photon field is a safe approximation for the quantum 
contribution in the late Universe. The purpose of the present work is to verify whether 
these quantum terms are compatible with the well-known classical cosmological solutions 
for the different epochs of the history of the Universe. The method which will be used in 
what follows is based on the derivation of metric perturbations over the given cosmological 
solution, in the theory with anomaly-induced effective action of gravity. 

The investigation of stability of the physically important solutions in a more general 
gravitational theories has a long and interesting history, starting from [15], where the sta- 
bility of Minkowski quantum vacuum has been explored for the first time (see also [T6]). 
Similar program for the theory of quantum gravity at non-zero temperature has been car- 
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ried out in [T7]. The stability of Minkowski has been further studied in [18]. The stabihty 
of de Sitter space for both semiclassical theory and quantum gravity, within different 
quantization schemes has been recently considered in [19]. Furthernmore, the stability of 
de Sitter space has been recently discussed in f20] by an unconventional method of power 
counting in the IR limit, with an apparently negative result concerning the validity of 
the whole semiclassical approximation (see further references therein). Compared to the 
methods of exploring stability used in the papers mentioned above, the anomaly- induced 
effective action has two important advantages: safety and simplicity. As we have already 
mentioned above, the anomaly-induced action is picks up the most important non-local 
part of effective action. As far as it is based on the conformal anomaly of massless fields, 
it is closely related to both UV and IR limits of the theory and hence is independent on 
whether we use usual in-out or more complicated in-in formalism. Also (see next sec- 
tion) the anomaly-induced action is very simple, it is given by a compact and explicit 
non-local expression which can be easily made local by introducing two auxiliary scalar 
fields. Historically, the anomaly- induced effective action was the theoretical basis of the 
first cosmological models with quantum corrections [21], Starobinsky inflation [7] and first 
derivation of cosmic perturbations in inflationary model [22]. As we will see below, the 
use of anomaly-induced action enables one to explore stability not only for Minkowski or 
de Sitter spaces, but for a wider set of classical cosmological solutions. 

In our investigation of stability we will restrict consideration by the gauge invariant 
part of metric perturbations related to the gravitational wave. The solutions of our interest 
include the radiation and matter - dominated Universes, the late Universe dominated by 
the cosmological constant and the stable phase of the modified Starobinsky inflation 
[231 El]. 

The paper is organized as follows. In Sect. 2 present a brief review of anomaly- 
induced action of gravity and of the corresponding cosmological solutions. Sect. 3 is 
devoted to the derivation of gravitational waves equation on an arbitrary cosmological 
(homogeneous and isotropic) background. In Sect. 4 we explore the stability of classical 
solutions by means of approximate analytic method and in Sect. 5 discuss the spectrum of 
metric perturbations. Sect. 6 is about the behavior of gravitational waves in the Modified 
Starobinsky model of inflation. Finally, in the last section we draw our conclusions. 

2 Effective action induced by anomaly 

The covariant form of anomaly-induced effective action of gravity [25l |26] is the most 
complete available form of the quantum corrections to the gravitational action in four 
space-time dimensions. The application to cosmology has been considered in [211 12Z] and 
led to the well-known Starobinsky model of inflation [7] (see also more detailed description 
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in [28] and consequent development in form of Modified Starobinsky model [23| 12^)- 
The anomalous trace of the energy momentum tensor is given by the expression [HI [2] 



< TH > 

fj. 



{wC^ + bE + cUR) 



(5) 



where the coefficients w, h and c depend on the number of active quantum fields of different 
spins, 
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No Ni/2 N, 
120 20 10 
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360 



lliVi/2 ^ 31 A^i 



360 



180 



(An) 



No iVi/2 

180 30 10 



(6) 
(7) 
(8) 



It is easy to see that these coefficients are nothing else but the /3-functions for the 
parameters 01,2,3 in the classical action of vacuum (jl]). Due to the decoupling phenomenon 
|14j . the number of active fields can vary from one epoch in the history of the Universe 
to another. As we have already mentioned in the Introduction, the present-day Universe 
corresponds to the particle content with A^q = ^1/2 = and Ni = 1. 

The anomaly-induced effective action Find represents an addition to the classical action 
of gravity, and can be found by solving the equation 



2 S Tind 



<Tll>= {uC^ + bE + cUR) . 



The covariant and generally non-local solution can be easily found in the form 
3c + 26 



36 



d!^x\/ —g{x) R^{x) 



(9) 



(10) 



+ j d^x^-g{x) j dS^/-9{y) {E -'^nR),G{x,y) 
where G{x, y) is a Green function for the operator 



^C^-^-iE-^-BR) 
4 8 ^ 3 ^ 



A4 = + 2 R^'V^V, - ^ i?n + i (V^i?) 



Finally, one can rewrite ( ITOl) in the local form by introducing two auxiliary fields and 
ip [29] (see also [30] for an alternative albeit equivalent scheme), 

3c + 26 



ind 



~ Sc[9^J.ly 

+ V 



36 



j d''x^-g{x)R\x) + j d*x^-g{x) |^^A4<^- 



w 



2 ^ 3 ^ 2v^ 



+ 



w 



'IV. 
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The expression (fTTi) is classically equivalent to (fTOl) . because if one uses the equations 
for the auxiliary fields and if), the nonlocal action (11 01) is restored. Consider now 
the background cosmological solution for the theory with the action including quantum 
corrections, 



total 



-Ml 



d'^x \/ —g -R + r 



121 



where Mj, 



l/lGvrG is the square of the Planck mass, and the quantum correction F is 
taken in the form ( ITTl) . 

Looking for the isotropic and homogeneous solution, the starting point is to choose 
the metric in the form g^^y = a?{ri)g^y, where rj is conformal time. It proves useful to 
introduce the notation a = In a. The theory includes the equations for the three fields, 
namely for Lp, ip, and a. For the sake of simplicity we will consider conformally flat 
background and therefore set g^u = rj^^. 

Equations for if and ip have especially simple form 



-9 



2 ^ 3 ^ 2v^ 







w 



2a/=6 



0. 



By using the transformation laws for the quantities in the last expression, one can obtain 



4 5 



(13) 



V^{E - 3 n^) = V^iE - -an + 4A,a) . (14) 

Taking into account our choice for the fiducial metric g^i, = rj^^, all the terms in the r.h.s. 
of the last equation are equal to zero and we arrive at the following equations 

□ V + Sttv^ □ V = , ^2^ = 0. (15) 

The solutions of ( lT5l) can be presented in the form 

</3 = -StTV^^ (T + (y9o , ijj = IpQ. (16) 

where □ is the fiat-space D'Alembertian and (po, iJjq are general solutions of the homoge- 
neous equations = , n"^ ipQ = 0. 

There is an obvious arbitrariness related to the choice of the initial conditions for 
the auxiliary fields ip, ip. However, replacing Eq. ( [T6|) back into the action and taking 
variation with respect to a we arrive at the unique equation for a. It proves useful to write 
this equation in terms of a(t) and the physical time t, with derivatives denoted by points. 
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Another useful variable is of course the Hubble parameter, H{t) = a{t)/a{t) = (T{t). 
Then, we obtaiE§| 

a 3 ad a ( 46\ da f 2b 

— + —1^ + — - 5 + — — -2A;1 + 



a a^ a^ \ c / a^ \ c / a^ 

Ml (d a 2A \ pO 



57rc \aa^ 6 j ca-^ 

In the r.h.s. of the last equation we have included the contribution of matter with constant 
Pm and explicitly shown dependence on a{t). As far as we deal with the trace of the 
generalized Einstein equations (linear combination of generalized Friedmann equations) , 
the contribution of radiation does not show up, but it can be easily restored if we switch 
to an equivalent (OO)-component [HI [121 IZ]- 

There are several relevant observations we have to make about the solutions of Eq. 
(fT7|l in different physical situations. First of all, in the theory without matter, when 
Pm = 0, there are two exact solutions, namely 

a{t) = ao ■ exp{Ht) (18) 

where |2l] 

1/2 



Mp I 647r6 A 



As far as the cosmological constant is quite small compared to the square of the Planck 
mass, A ^ M|,, we meet two very different values of H (here A > 0) 

It is easy to see that the first solution with He is the one of the theory without quan- 
tum corrections, while the second value Hs corresponds to the inflationary solution of 
Starobinsky [7]. The sign of b is always negative, independent on the particle content, see 
Eq. ([7]). Let us remark that the particle content which we deal with here, Nq, N1/2 and 
Ni , corresponds to the degrees of freedom contributing to the vacuum effective action 
and has nothing to do with the real matter content of the universe. 

Second, the stability properties of the solutions fl20|) depend on the sign of the coef- 
ficient c, that is on the coefficient of the local J y/—gR^-teTm [71 123] . The inflationary 
solution Hs is stable for a positive c and is unstable for c < 0. The stability of the low- 
energy solution He requires opposite sign relations for c. It was shown in [24] that the He 



^We have included here k = 0,±1 for the sake of generahty, but the rest of the paper will be only 
about the k — case. 
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solution with c < is stable with respect to the small variations of the Hubble parameter 
(or, equivalently, of cr(t)). In the original Starobinsky model of inflation [7] the particles 
content corresponds to the unstable inflation and the initial data is chosen such that the 
Universe is asymptotically approaching the radiation-dominated FRW solution. In order 
to better understand the situation, let us replace the corresponding FRW solution, e.g., 
a{t) ~ t^/^, into the equation ( ITTl) . It is easy to see that the classical part, composed by 
Einstein and matter terms, do behave like while the quantum corrections, which are 
given by higher derivative terms, behave like This means that in the unstable phase 
the quantum terms do decay rapidly, such that the classical solution a{t) ~ t^/^ is an 
excellent approximation to the solution of Eq. (fTTll in the corresponding epoch. It is an 
easy exercise to check that the same is true for the radiation-dominated a(t) ~ t^^^ and 
cosmological constant-dominated epochs too. An alternative form of these considerations 
can be found in [31]. Indeed, all the arguments presented above are valid for the dynamics 
of the conformal factor only. The next sections will be devoted to the stability of the same 
classical solutions with respect to the tensor metric perturbations (which we identify, for 
brevity, as gravitational waves). In the theory with quantum terms such as (jl]) and ( ITTi) . 
the equations for the gravitational wave have fourth derivatives and hence they represent 
a real danger for stability of the classical solutions. 

In the modified version of Starobinsky inflation [231 El] the Universe starts with the 
stable inflation, that can be provided by choosing the supersymmetric particle content 
and Ni . Then the exponential inflation slows down due to the quantum effects 
of matter fields (mainly s-particles) and at some point the heavy s-particles decouple from 
gravity and then the Universe starts the unstable inflationary phase. The advantage of 
this version is that it does not depend on the choice of initial conditions. The stability of 
the stable inflation with respect to the metric perturbations will be explored in Sect. 6. 

3 Derivation of the gravitational waves 

Let us derive the equation for the tensor modes of metric perturbations. First we rewrite 
the action in a more appropriate way and then derive linear perturbations for the tensor 
mode. 

3.1 Total action with quantum terms 

It proves useful to present the action (fT2|) of a more useful way. After performing some 
integrations by parts, it can be cast into the form 




(21) 
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with 

5 

L = Y,fsLs (22) 

s = 

= [/o i? + /i R^^^'Ro.p^.u + /2 R'^Ro.^ + hR^ + /4 v^R + h ^^v] , 

where the /-terms are defined as 

/o 
/i 



_ 

IGvr ' 



h + uj uj 



2V-6 



2V-6 



/2 

/s 

/4 
/5 



UJ + 2h u; , 

-2ai - 4a2 H -j= tp; 

ai 3c + 26 36 + a; u , 

h ao := H ^= "0; 

3 36 6v^ 6v^ 



(23) 
(24) 

(25) 

(26) 

(27) 
(28) 



where the coefficients ai^2 are the same defined in 



3.2 Perturbation equations 

Using the conditions (with /i = 0, 1, 2, 3 and i = 1, 2, 3), 

di h^^ = and hkk = , 



(29) 



together with the synchronous coordinate condition /i^o = 0, we introduce metric per- 
turbation in the equation (122|) as follows 



(30) 



Here (7°^ = {1, —5ij a?{t)} are the background cosmological solutions. In this way one 
can arrive at the following expressions for the bilinear parts of the partial Lagrangians 
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fromEq. ([22]): 
Lo - 



a^fo 



( 2H^ - 2H 



hh + AHhh + -h^ - -^r- 



hh(8HH + 1QH^ 
H'h-2Hh'^^^ 



hh[AH^ 



+ 2h^- + 



a'f2 



+ AH 
+ h'^ + AHhh + 

\ h^ 

hh{ 12HH + 2AH^] {5H+^H 



+ 0{h^ 



h^ 3H^ + QHH^ + - 



- h^ 3H^ + 9HH^ + 9H^] -hh{AH + 6H^] + — + -Hhh + 



l/V2/i\2 1,, 

--{h + 3Hh 



4 



Qa^ h[H + 2H' 



Hh - 3H^h 
3H^ 



V^h 



+ 2hh + 



4 

+ 0{h^ 



+ 8Hhh + -h^ — 



^(H + 2H^ 



+ 



^{^ + 3H^]h' 



+ 0{h'), 

^ + 3H^ h^ ^[(p^ 3eS) hh + (iSH'^Lp + AHip^ hh + 
1 



Lfh"^ — —H(^h'^ — Lfhh 



3Hip)—hV^h +0(h^] 



+ 0^ 



4 



-Hh' 
4 



lHhh--hh~-Xv'h 
3 3 6a2 . 



+ 0{h^ 



(31) 



One can perform a comparison of these equations with the ones known from the hterature. 
A very similar expansion was obtained by Gasperini in [32] in order to explore metric 
perturbations in the pre-Big-Bang inflationary scenario. Despite the physical motivations 
of the pre-Big-Bang inflation are quite different from our case, when we intend to consider 
the semiclassical gravity with the action induced by anomaly, the formulas are mainly 
equivalent and we find a perfect correspondence between our results and the ones of p|2]. 
On the other hand, we can successfully compare part of the expressions presented above 
with our own previous calculation of the same equations for the H = const case in [T3] . 

In order to get the equation for linearized tensor perturbations, one has to omit all 
higher order terms 0{h^) in the expressions ( 1311) and then proceed by taking variational 
derivative with respect to h^^^. The next step is to use the solutions (|T6]) for the auxiliary 
fields. We fix the ambiguity in these solutions by choosing the simplest zero option for the 
conformal terms which are not controlled by conformal anomaly and set ipQ = ipo = 0. 
As a consequence we have 



in 



hH ^ 



57r 



hH . 



hH and ip= -SnV-b H . (32) 
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These relations for the background must be replaced into the equations for tensor pertur- 
bations, but we prefer to keep (/^-dependent form, for the sake of generality. After all, the 
equation for tensor mode can be cast into the form 



+ 



+ 



2/1 + y) h + 
9 X 



3iJ 4/1 + /2 +4/1 + /2 



h + 



3if^(6/i + ^-4/3 



+ H I6/1 + /2 + 6iJ /i - /s + 2/1 + - [h + fo + U^) + -UH^ - -U^' 



4/1 + h 



H[Ah-Qh]-2lHH[-h + 2h)-H[-h + Qf, 

3 •„■ 



3if2 4/1 + -/2 - 4/3 - ^H\h + 4/3) +HUh + -h) + -fM^H^ + H 



H{^h(p + \h - 2/50') + \h -\h^v\ h - [H{Ah + /s) + 4/1 + h] 
5UH +U V - {36HH^ + 18H^ + 24HH + 4 # ) (/i + /a + 3/; 



- M(32/i + 36/2 + 120/3] - 8//(/i + /2 + 3/3] - ff' (4/1 + 6/2 + 24/3 

- 4H(h + f2 + 3/3) - 9U^(h^ + + /4¥3f3ff2 + 5h) - H^fsh + I2/2 + 48/- 



6' 



2H + 



1 • 



/7^(4/i+4/2 + 12/3) +/J(2/i + -/2) +2if(/i + /2 + 3/3 



2/1 + 2/2 



0, 



(33) 



Finally, if we take H = constant in f l33|) . we find the result which perfectly fits the 
one of IT3l. 



4 Stability analysis 

Now we can start to deal with our main task and see whether Eq. ( 133|) indicated that 
there is a stability of the cosmological solutions or not. As a starting point we remember 
that a linear dynamical system with constant coefficients is stable when all poles, i.e., all 
roots of its characteristic equation have negative real part, i.e., they are on the left half 
complex plane. We can analyze this problem in two ways: numerically and analytically. 
We will start with an approximate analytical analysis. 
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4.1 Semi- analytical analysis 



Let us start by rewriting the terms in (1331) by using the plane wave representation in fiat 
space section, 



h 



-n 



—n 



aity ' 
h 



h 



-n 



a{tY 
h 



n 



Then the equation for tensor perturbations can be presented as follows 

h'h +hh + h2 h + hi h + ho h= 
where we used the notations 

3iJ(4/i + /2) +4/1 + /2, 



2A + f 



b2 



n 



f' 



4/1 + /2 K + 3if^ 6/1 + ^ - 4/3 + ^ I6/1 + 7t/2 )+6H(^-f: 
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+ 2 A + ^ (/2 + /O + U^) + \fAH^ - ^/50' , 



hi = [if(4/i + /2) +4/1 + /2] ^ + if(4/i-6/3) -21m(^/2 + 2/3 

- Hi^-h + 6/3) + (4/1 + \f2 - 4/3) - 9iJ3(/2 + 4/3) + (4/1 + 
+ y,^[m^ + h)+ H{3U^ + - 2/50') + V -^/5W , 

60 = 5UH V? +/4 'V^' - {36HH^ + 18H^ + 2AHH + A H ) (^fi + f2 + 3/3) 

- m(32/i + 36/2 + 120/3) - 8if (/i + /2 + 3/3) - (aJ, + 6/2 + 24/; 

- AH (h + f2 + 3h)- 9/40 (h^ + HH] + /40 f 3if2 + hH 
r3 



(34) 
(35) 

(36) 
(37) 

(38) 



(39) 



H' ( 8/1 + 12/2 + 48/3 ) + hv' {^H^ + \h) + \hmv - \h(p^ + V 



+ /o 
1 



2i7 + 3/7' 



if^(4/i+4/2 + 12/3) +//(2/i + -/2) +2if(/i + /2 + 3/3 



+ 77 (/2 + /40 + /O + 3/4^/0) --/50' 



J 



1 



+ 2H 



n 



(40) 



2 •' 3' 

The auxiliary field and its derivatives should be replaced according to Eq. 

So, we have to analyze the equation ( 135|) . with the coefficients hk, {k = 0,...,3) 
given by equations (137|) . ( 138|) . ( |39|) and (140|) . One can easily reduce this fourth-order 
equation to a system of four first-order equations. Making the new change of variables 
we introduce 



ho = h , hi = ho 



h , /12 



h. 



h, /i3 



h. 



h . 



(41) 
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Rewriting the differential equation, we arrive at 

&3 + &2 h2 + bihi + bo ho), 



CO 


1 


~h 


h2 = 


hs, 




h2 , 


ho = 


hi. 



We can rewrite the hnear system of four equations given above in a matrix form to 
compute its eigenvalues and eigenvectors. Thus, we can write in simplified form. 



where k 



hk ^ A[hi, 
0, 1, 2, 3 and the matrix A 



(42) 



A[. has the form 








1 


















1 


















1 




v 


do 


di 


d2 


ds 


1 



A 



Here we called dk = ~bk/hi. 

If the coefficients of the matrix A would be constants, the problem of stability could 
be solved immediately by deriving the eigenvalues of A. However, the same method is 
know to work even for the non-constant matrix A. The reason is that we are looking 
for the asymptotic stability related to the exponential time behavior of hk- As far as the 
coefficients of the matrix A have weaker time dependence, we can neglect this dependence 
and treat ^ as a constant matrix. It is easy to see that this condition is satisfied for 
radiation and matter - dominated backgrounds and can be used also for the cosmological 
constant - dominated background, because A is actually very small. In other words, our 
constant-A approximation means that we are looking for the dynamics of perturbations 
which is stronger than the expansion of the Universe in a given epoch. As far as we 
are interested in the consistency of classical cosmological solutions and in avoidance of 
dangerous run-away type solutions, this is definitely a very reliable approximation. 

So, the next task is to find the eigenvalues of A and hence we consider 



det 



The algebraic equation 



1 


1 








\ 





-A 


1 













-A 


1 






di 


^2 




A) / 


-d^X" 


— d2 


A^- 


diX^ 


- do 



(43) 



(44) 
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After some algebra (see Appendix for details), we can reduce the above equation to the 
following form 

+ ^1 ^ + 6 = . (45) 

The most important quantity is 

3 



6^^ /6 



+ ■3 



(46) 



The value of A, obtained by using the Cardano formula, and all notations used here, are 
explained in Appendix. Eq. fl75]) . will tell us the nature of these roots. We can distinguish 
the following distinct cases: 

1. A < 0: Then the three roots are real and distinct and can be, 

• All negative roots: stable. 

• Some positive root: unstable and instability generally increases with increasing 
number of positive roots, in a sense one needs more severe initial conditions to 
avoid instability. 

2. A = 0: the roots are real, and two or three are equal. Then, 

• All negative roots or with negative real parts: stable. 

• Some root with a positive real part: unstable and this instability increase with 
increasing number of such positive roots. 

3. A > 0: one real root and two complex roots, 

• All negative roots or with negative real parts: stable. 

• Some root with a positive real part: unstable and this instability increase with 
increasing number of positive roots. 

In the case of equation (|33|) one meets the following values for (j46l) . 
6 = — + /3 and 6= h7 + ^ 



a 



5 1 / \ 

- p; 1 = -{q'^ — 4p^ + Apr) and /3 = — r , (47) 
2 8 V / 



p = di + do] q = — h dl and r = ^ H ^ h an • 



Let US remember that 6a,./&4 = ~dk, where bk are given by (137|) . (!38|) . (!39|) . ( |40|) . Hence 
one can expect that the expression for A from the equation (146|) will be rather complex, 
requiring numerical analysis. 
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By performing such an analysis for the three cases of our interest, namely for expo- 
nential expansion, radiation and matter epochs, we find 

1. Exponential expansion. When we choose ai < Owe found A < 0. This is consistent 
because when we analyze equation (jUj) directly, we find all eigenvalues to be real 
and negative. So, we have the stability in this case, exactly as we could expect from 
comparison to the inflationary case |24j . 

However, if we choose ai > 0, then we find A < too. But analyzing equation 
(jS]) directly by numerical method (this means deriving the roots numerically by use 
of Mathematica software [33]), we find three negative and one positive eigenvalue. 
So, we can observe the instability in this case. 

Let us remark that the sign of Oi defines whether the massless tensor mode in the 
classical theory is a graviton or a ghost [5] . From this perspective our result means 
that the stability property of the theory with higher derivative classical term 
and quantum correction (ITOl) is completely defined by classical part (I3l) and, quite 
unexpectedly, does not depend on the quantum term ffTOl) . Is it a general feature or 
just a peculiarity of the de Sitter background solution? Let us consider other cases 
to figure this out. 

2. Radiation. When we choose oi < we found A > 0. This is consistent because 
when analyzing the Eq. f H4|) directly, we find two real eigenvalues, which are both 
negative and also two complex eigenvalues with negative real parts. So, we have 
stability in this case. But if we take ai > it turns out that A < 0. Analyzing Eq. 
( l44l) directly, we find two negative eigenvalues and two positive ones. So, we have 
instability in this case. Again, the stability of the classical solution is completely 
dependent on the classical term (|3]). 

3. Matter. With ai < we find A > 0. This is consistent the direct numerical 
analysis of Eq. (jHj), because in this way we find two real negative eigenvalues and 
also two complex eigenvalues with negative real parts. So, we have the stability for 
ai < 0. However, if we choose ai > 0, we find A > 0, indicating instability. By 
analyzing Eq. ( I44l) directly one confirms this result, for we meet two real eigenvalues 
(one negative and other positive) and two complex ones, both with negative real 
parts. 

As a result of our consideration we can conclude that there is a stability for Eq. ( l33|) . 
if and only if ai is negative. Taking into account the mentioned feature of classical higher 
derivative gravity, we see that the linear (in) stability of tensor mode in the classical higher 
derivative theory ([3]) completely defines a linear (in) stability in the theory with quantum 
correction (fTOj) . The quahtative explanation for this output is quite clear. The quantum 
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terms (fTO!) consist of two types of terms. The simplest one is the local i?^-term, which 
contributes to the propagator of gravitational perturbations on flat background, but not to 
the one of the tensor mode (see, e.g., [3] for detailed explanations and original references). 
The more complicated non-local terms are at least third order in curvature, and hence do 
not contribute at all to the propagator of gravitational perturbations on flat background. 
Indeed, we are interested in the perturbations on curved cosmological background and not 
on the fiat one. However, the typical length scale related to the expansion of the Universe 
is defined by the Hubble radius and are much greater than the length scale of the linear 
perturbations we are interested in here. Therefore, the stability of theory under such 
perturbations, in a given approximation, is the same as for the fiat background and for 
the algebraic reasons explained above, there is no essential role of the anomaly-induced 
quantum terms f lTOj) here. 

4.2 Numerical analysis 

In order to ensure that our qualitative and analytic consideration of the stability is correct, 
let us present the analysis of the stability of the differential equation fl33|) by means of 
numerical methods, using the software Mathematica |33j . 

We tested the same three relevant cosmological solutions, namely exponential expan- 
sion, radiation and matter. In all cases the initial conditions of quantum origin were 
taken, for the sake of simplicity. For the mentioned three cases one can meet both stable 
and unstable solutions, as they are shown in Figure 1 and Figure 2. One can easily see 
that the solutions where it was adopted oi < are always stable, as shown in Figure 
1. At the same time the solutions where it was adopted an opposite sign, ai > 0, are 
unstable, as shown in Figure 2. 



Figure 1: Graph of perturbation h{t) in function of time analyzed in the cases for a{t) = 
Ooe^"*, a{t) = aot^^'^ and a{t) = aot^^^ respectively, with the initial conditions = 
h = = h= where we adopt ai < 0. Stable behavior. 
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Figure 2: Graph of perturbation h{t) in function of time analyzed in the cases for a{t) = 
aoe^°*, a{t) = aot^^"^ and a{t) = aot^^^ respectively, with the initial conditions = 
h = s/^^h = h= where we adopt ai > 0. Unstable behavior. 



5 The spectrum of gravitational waves 

Before starting to work with the differential equation (133|1 . let us remember some standard 
notions concerning the spectrum of gravitational waves. For the sake of simplicity we start 
from the inflationary background case. 

Consider the case of gravitational wave on the background of inflationary solution in 
the theory without quantum corrections. We have 

■■ ah ( 2a^ , ^ 

h + — \h = 0, (48) 

where we assume a{t) = a^e^^* and also set oq = 1. The vacuum state is well known 
for this theory (see, e.g., |2j), being 

Mx, r^) = Mr/) e±'"■^ /,(r^) = ^ , (49) 

v2n 



where n is called the wave number vector. In order to study the dynamics of h{t, x), it is 
necessary to make a Fourier transform, 

Kit) = I h{t,^)e^--d'x. (50) 

We will need the total square of the amplitude, namely 

h\t) = J hl{t)d?n. (51) 

The above equation can be rewritten in the form 

h^{t) = 47r y hl{t)n^dn = ^tt ^ hl{t)n^d\nn = j Pn{'t)d\nn, (52) 

where 

Pnit) = hl{t)n\ (53) 
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The last quantity is called the "square of the power spectrum". It shows us how the 
amplitude of gravitational waves vary in a range of In n to In {n + dn). 

Starting from this point one can find the power spectrum for theory with Eq. ( HHl) 
for the gravitational wave and then apply in our model with Eq. fl33l) . For this end one 
has to square the value of the gravitational perturbation at a given time and for a given 
wave number n. We will vary this n for a fixed t and simultaneously solve our fourth- 
order differential equation numerically. After this, we linearize the graph by plotting the 
relation 

In n^hl{t) X In n. (54) 

As a result we obtain the linear proportionality coefficient, which will be denoted as k and 
called spectral index. Then we have, -P^(t) oc , i.e., it is proportional to the spectral 
index. It is the power spectrum that will tell us how the amplitude of the perturbations 
depends on the wavelength. 

Let us adopt the initial conditions of quantum origin, 

hooc-^, hoc^-, hoc^, hoc-^. (55) 



Then, performing the above procedure for the case of inflation, represented by the equation 
fj48ll . we found, 

k ^ 0.01 . (56) 



Now, applying the same procedure for our model of equation ( 1331) . we find, in the case of 
exponential expansion without quantum corrections (this means ai = a2 = 0), 

k = 0.00894167 ^ 0.01 (57) 

where we use the units with the Planck mass equal to 1 in the equation ( 133|) and the Eqs. 
fl28|) . for perform the numerical analysis. This value is very close to the one which can be 
found for infiation when we calculate analytically, it should be 0, as we can see on 
( l56l) . Thus the model we are considering returns to the well known infiation case when 
we take ai = 02 = 0. The fiat or almost fiat spectrum which is found occurs because the 
stable version of the anomaly-induced infiation develops de Sitter phase with H = const. 

Now we are able to analyze the general quantum case with ai 7^ 0. Remember that 
the action for the vacuum is given by ([1]), where ai is most important for the evolution 
of tensor perturbations. Below we consider three separate cases and present the results 
of numerical analysis. 
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5.1 Inflation, or exponential expansion 

The variation of the power spectrum with respect to the change of ai (positive and 
negative) is presented in Table 1. 
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Table 1: Getting k through the variation of ai > in case of a{t) = UqC 
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Table 2: Getting k through the variation of ai < in case of a{t) = aoe 



These two tables show that when the value of ai is negative and is decreasing, the 
values of k also decrease while being negative. That is, as we increase ai, the amplitude 
of gravitational waves decrease. In the case of positive oi an opposite situation happens, 
namely when we increase its value, the value of k is also increasing. This happens because 
ai > does not represent a stable solution of ( 133|) . 

These results may be compared with the computation of the spectral index for grav- 
itational waves in the context of the pre-big bang scenario, based on the string effective 
action at tree level [34j. In this case, the spectral index is positive with the increasing 
spectrum, while the de Sitter inflation predicts a fiat spectrum and power law inflation a 
negative spectral index, that means the decreasing spectrum. 

Now, using the data obtained in [35], [36], [37] and using the software CMBEASY 
[38] [39] and [10] we can obtain the graph of the spectrum of anisotropy of the cosmic 
microwave background (CMB) due to gravitational waves for the cases specified above. 
In these graphs there are two essential quantities, namely the spectral index k and the 
content of matter in the universe. The plots for the exponential expansion are shown in 
Figure 3. In all graphs we use 02 = 0. However, k value is the same for any 02 value. 
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Figure 3: Behavior of the anisotropy spectrum for the three cases with k = 0.01. These 
three solutions have the same to ai = 0. That is, they are degenerate. 

5.2 Radiation 

By performing the same procedure of solving numerically the differential equation (!33|) 
and the same linearization as used before, we can test other cases, like radiation where 
we have a{t) = a^t^^'^ and matter, where a{t) = aot'^^^. 
For radiation, we have the results presented in Table 3. 
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Table 3: Getting k through the variation of ai > in case of a{t) = aot^^'^. 
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Table 4: Getting k through the variation of ai < in case of a(t) = aot^^'^. 



With these values we obtain the graphs shown in Figures 3 and 4, for the anisotropy 
spectrum, where we take the value ai = (Figure 3) and the extreme value ai = 1 
(Figure 4) from the tables presented above. 
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5.3 Matter 



Finally, we analyze the case for matter-dominated epoch, with a{t) — aot^^^ , where we 
get the values presented in Tables 5 and 6. 
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Table 5: Getting k through the variation of ai > in case of a{t) = aot^/^. 
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Table 6: Getting k through the variation of Oi < in case of a{t) = aot^/^. 



With these values we obtain the graphs for the anisotropy spectrum presented in 
Figures 3 (ai = 0) and 4 (ai = 1). 

We can observe that the graphs of the spectrum of CMB anisotropy are nearly identical 
in the case of radiation and matter (Figure 3), but both are different from the case of 
inflation. The cases of matter and radiation in Figure 4 are degenerate. 
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Figure 4: Behavior of the anisotropy spectrum for gravitational waves in de Sitter (black 
line with k = 0), inflation (red line with k = 0.79) and Radiation/Matter (blue line with 
k = 0.3 for both). All this is for ai = 1. 



6 Quantum effects of massive fields and tempered 
infiation 

It was shown in [23| that taking the quantum effects of massive fields into account leads 
to the tempered form of inflationary solution Hs in ( 120|) . The solution which emerges 
after we derive the corresponding effective equations can be very well approximated by 
the formula 

a(t) = e'^W , a{t) = H,t - ^ ft' , H{t) = Ho-^ft, (58) 

where / is a small dimensionless parameter, which is at least as small as 10~^. The plot 
of this parabolic a{t) this given in Figure 5. This stable phase of inflation is supposed to 
last until H{t) becomes comparable to the energy scale of supersymmetry breaking and 
the transition to the unstable phase of infiation occurs. 
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Figure 5: Evolution of cr{t). 

We can use the solution (155]) in our differential equation (155]) for the tensor perturba- 
tion. Solving this equation numerically and doing linearization as it was explained above, 
we arrive at the Tables 7 and 8 for the variation of k with respect to ai in the case of 
Modified Starobinsky model of inflation. 
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Table 7: Getting k through the variation of ai > in case of Eq. (158]) . 
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Table 8: Getting k through the variation of ai < in case of Eq. (158]) . 



We can see that the table is almost identical to the one for the case of usual inflation, 
described in Table 1. This is because the parameter / is small 10~^ and the desacceler- 
ation of inflation happens slower compared to the typical scale of tensor perturbations. 
If we consider this parameter being increased, the two results begin to diverge from each 
other. The values of k increase as the value of / increase too. 

7 Conclusions 

The classical action of gravity which is necessary to guarantee the renormalizable theory 
of matter fields on classical metric background goes beyond the conventional Einstein- 
Hilbert term, for it includes also the cosmological constant term and higher derivative 
terms. The presence of higher derivatives is usually associated to the problem of higher 
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derivative unphysical ghosts, but this issue does not represent a real problem for the theory 
when gravity is only an external background and there is no unitarity condition for the 
gravitational 5*- matrix. In the case of external gravity the necessary consistency conditions 
should be the existence of physically acceptable solutions and their robustness. The last 
criteria means that these solutions must be stable, at least, with respect to the small 
perturbations of the metric variables. In case of gravity the higher derivative terms can 
be especially dangerous for the stability of classical cosmological solutions, such as matter 
and radiation - dominated Universes and the Universe dominated by the cosmological 
constant term. The results found here shows that the detection of the gravitational wave 
contribution to the spectrum of anisotropy of the cosmic microwave background radiation 
(which may be possible with the Planck satellite) could lead to important bounds on the 
possible quantum effects analysed in this paper. 

On the top of the classical vacuum terms there are also quantum corrections to it. 
One of the best available forms of such corrections is the anomaly- induced effective action 
of vacuum, which is especially efficient as an approximation to the present-day Universe, 
when all quantum fields except the photon have already decoupled from gravity and the 
quantum contribution of photon fits perfectly to the anomaly-induced approach. 

We have used the anomaly-induced action of gravity to explore the behavior of ten- 
sor perturbations and especially to verify the stability of classical cosmological solutions 
with respect to such perturbations. The consideration was based on explicit derivation 
of gravitational wave equations in the theory with anomaly-induced quantum corrections 
and on the use of both analytical and numerical methods to perform the detailed analysis 
of these equations. The main conclusion of our work is that the stability conditions are 
essentially related to the sign of the Weyl-squared term in the classical action of vacuum 
and do not manifest any essential dependence on the quantum contributions. The quali- 
tative explanation of this result is that the anomaly-induced action has a structure which 
prevents it to give essential contributions to the linear equation for tensor perturbations. 
Let us note that the situation may be very different for the density perturbations, where 
the sign of the i?^-term may be most relevant and for the non-linear perturbations, where 
the mentioned qualitative arguments simply do not hold. 
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A Roots of the quartic polynomial equation 

Here we present the analysis of the fourth-order polynomial which was used in Section 4. 
Consider the polynomial, 

+ + rfaA^ + diX^ + do = . (59) 

Let us follow the classical Cardano method. To eliminate the cubic term make the sub- 
stitution A = X — ^. We get 

x'^ + px^ + qx + r = , (60) 

where 

P = i'i + d^^ (61) 

. = + (62) 

3c?q dod'o dodi , , , 

- - +^-^ + do. (63) 



Now we rewrite fl^U]) as 



256 16 



x"^ + px^ = —qx — r 



by completing the square and adding the y term. We obtain 

y^ + ay^ + /3y + j = 0, (64) 

where 

a = —p, (65) 
(3 = 2p2-r, (66) 

7 = -{q^ - Ap^ + Apr) . (67) 
8 



In order to eliminate the second degree term in (l64l) . make one more substitution y = 
z + m, where m = —a/2. 

^' + ^1^ + 6 = 0, (68) 

where 

6 = + (69) 

6 - ^ + ^^- (70) 
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For the equation (168!) one can already apply the Cardano formula and find 



A 



The solution is 



6 

2 



or 



2 2 V 2 2 

As far as we know the values of z, m, a and p, it is possible to find the roots A. 



(71) 



(72) 



(73) 
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